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We analyze the change in free energy, internal energy and entropy due to the presence of a heavy 
quark anti-quark pair in a QCD heat bath. The internal energies and entropies are introduced as 
intensive observables calculated through thermal derivatives of the QCD partition function contain- 
ing additional static color sources. The quark anti-quark internal energy and, in particular, the 
entropy clearly signal the phase change from quark confinement below and deconfinement above the 
transition and both observables are introduced such that they survive the continuum limit. The 
intermediate and large distance behavior of the energies reflect string breaking and color screening 
phenomena. This is used to characterize the energies which are needed to dissolve heavy quarkonium 
states in a thermal medium. Our discussion supports recent findings which suggest that parts of 
the quarkonium systems may survive the phase transition and dissolve only at higher temperatures. 
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I. INTRODUCTION 



This is the second part of our discussion of thermal 
modifications of the strong forces in finite temperature 
QCD (for a detailed introduction to this subject and 
further references see 0, 0)- At finite temperature, 
T ^ 0, the free energy of a static quark anti-quark pair 
HL H- separated by distance r, is an important tool to 
analyze the in-medium modification of the QCD forces. 
Similar to the free energies also the internal energies have 
recently been introduced 0, |(| and are expected to play 
an important role in the discussion of quarkonia binding 
properties 0, 0, S EH E3, E3, E3, EM E3l More- 
over, the structure of these energies at large distances and 
high temperatures is important for our understanding of 
the bulk properties of the QCD plasma phase, e.g. the 
screening property of the quark gluon plasma |lcill9| |. the 
equation of state [2(J, |2l| . They also provide important 
input to the construction of effective models based on 
properties of the Pol^akov loop HEplllllllI^. Up 
to quite recently [H Ell I2I IM Wi l28|"most of these dis- 
cussions concerned the quark anti-quark free energies in 
quenched QCD. Several qualitative differences, however, 
are to be expected when changing from free to internal 
energies and/or when taking into account the influence 
of dynamical fermions. The difference between free and 
internal energy arises from non-trivial entropy contribu- 
tions [HQ- Moreover, in QCD with light dynamic quarks 
the large distance behavior of the strong interaction will 
show a qualitative different behavior due to the possibil- 
ity of string breaking. 



Properties of the finite temperature quark anti-quark 
free energies and the heavy quark potential at T = 0, 
U(r), have been discussed for 2-flavor QCD in Refs. 0, 
EH Ea , I27I |2S| . Recently some results have also been re- 
ported for 3-flavor QCD j2(|. Here we will continue our 
analysis of the fundamental forces of QCD at finite tem- 
perature. We analyze the partition function of 2-flavor 
QCD in the presence of heavy quarks and extract the 
quark anti-quark internal energies and entropies. This 
paper is organized as follows: In section^ we introduce 
the change in internal energy and entropy due to the 
presence of static quarks and anti-quarks in a thermal 
heat bath. We discuss their temperature dependence at 
large distances, in particular, their behavior in the vicin- 
ity of the transition, in section IIIII We finally discuss 
the qualitative and quantitative differences between free 
and internal energies in section lTVl and discuss their bind- 
ing properties with respect to quarkonium binding in the 
vicinity of the transition. Section Ivl contains our conclu- 
sions. Details on our simulation parameters, the lattice 
actions used in our calculations as well as details on the 
analysis of the quark anti-quark free energies are given 
in Ref. Q. 



II. PARTITION FUNCTION IN THE 
PRESENCE OF HEAVY QUARKS 

A. Free energy, internal energy and entropy 

As we are interested in the lattice formulation of QCD 
at finite temperature in thermal equilibrium, we consider 
the (Euclidean) path integral, i.e. we consider the parti- 
tion function of the QCD heat bath, 
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Z(T, V) = / dAd^d^>e 



-5[A,*,*] _ e -F(T,V)/T 



,(1) 



2 



where T (V) denotes the temperature (volume) and 
S[A, 'J, ty] the QCD action. We also investigate the corre- 
sponding system containing additional heavy quarks Q, 



functions the singlet (1), averaged (qq) and octet (8) color 
representations of the operator Or [W, W^] can be spec- 



Z (r, T,V) = / dAdVdV O r c) [W, W^] e 



-S[A,*,*] 



= e --Fo(r,T,V)/T 



(2) 



where 0^ [W, W^] denotes an operator which introduces 
static color sources representing quarks and anti-quarks 
separated by distances r = {r^} in a specific color 
representation c. A static color source appearing in 
[W, located at x is described by the thermal Wil- 
son line, 



W(x) = T exp 



-l/T 



dx A • A (x ,x) 



(3) 



the 



The expectation value of Or [W, W^}, i.e. 
point Polyakov loop correlation function (O r c ^ [W, W 1 )), 
is related to the change in free energy, Fq (r, T) = 
F(r, T, V) — F(T, V), due to the presence of static quark 
anti-quark sources in the heat bath, 

F (r,T) = -T\n(OW[W,W*])+TCa (4) 
= -TQnZo(r,T,V)-hiZ(T,V))+TCo , 

where Co can be fixed through renormalization |2^| . For 
instance, in the case of 2-point Polyakov loop correlation 



OP\W,W*] = -TrW(0) W\\r\) , (5) 

= -TrW(0)TrW^{\r\) , (6) 
9 

OW[W,W f ] = -TrW(0)TrW^{\r\) 
8 

-lTrW(0)Wt ( | r |). (7) 



The color singlet, averaged and octet quark anti-quark 
free energies, i.e. Fi(r,T), Fq q (r,T) and F s (r,T), re- 
spectively, have already been discussed extensively in 
quenched and full QCD HUHHIMIillillilS 

For the purpose of discussing internal energies 
(Uo(r,T)) and entropies (So(r,T)), we follow the con- 
ceptual approach suggested in Refs. @, [|| and consider 
thermal derivatives of the QCD partition functions intro- 
duced above, i.e. 



Uo(r,T) 



-T 



2 d F (r,T)/T 
dT 



(8) 



which leads to 



J 



U a (r,T) = -T 



[Or C) [W,W^]) 
Uo(r,T, V)-U(T, V) 



dT 



dT 



(9) 



r 



(c) 

Note here that the derivative of the operator O r with 
respect to temperature, dOr /dT, vanishes due to 10 . A 
similar relation can be derived for the entropies starting 
from 



So(r,T) 



dFg{r,T) 
dT 

S (r,T,V) - S(T, V) 



dTCp 
dT 



,(10) 



i.e. the observable TSo(r,T) could be calculated from 
the difference, TS {r,T) = U (r,T) ~ F (r,T), with 
U {r,T) and F {r,T) given in © and QJ. Wc navc 
also specified constant contributions, Co, which result 
from divergent contributions to the free energies and will 
control the internal energies and entropies at large quark 
anti-quark separations. Once the free energies are fixed 



through renormalization also the constant contributions 
to the internal energies and entropies are properly deter- 
mined. 

We note that Eq. @, and similar 1)11) |). open the pos- 
sibility for a direct calculations of Uo(r, T) and So(r, T) 
and define properly the quantities we aim to discuss here, 
i.e. the change in internal energies and entropies due to 
the presence of static quarks and anti-quarks in the QCD 
heat bath. Quite similar to the change in free energies, 
F (r, T) = F (r, T) - F(T), also the changes in internal 
energies and entropies, U {r, T) = U {r, T, V) - U(T, V) 
and S {r,T) = S (r,T,V) - S(T,V), are expected 
to behave like intensive observables and, in particular, 
will show no volume dependence in the thermodynamic 
limit. It, however, can no longer be assumed that the r- 
dependence of the quark anti-quark free energies are 



3 



given by the r-dependences of the internal energies © 
alone, i.e. we expect 

F a (r,T) = U (r,T)-TS (r,T) . (11) 

This indicates a quite complicated relation between free 
energies, internal energies and entropies and stresses the 
important role of the entropy contribution which is still 
present in free energies. In the limit of vanishing tem- 
perature, T — > 0, however, TSo(r,T — > 0) will vanish 
and 

F o {r,T = 0) = U o (r,T = 0) = V a (r) , (12) 

i.e. F D {r,T -> 0) and U {r,T -> 0) could be used 1 to 
define the corresponding energies at T = 0. 

Unfortunately, the non-perturbative formulation of en- 
ergies given in is still complicated and, in particular, 
it is complicated to be give n in a form suitable for lattice 
simulations |2lt 139 . l36l l37| . In the thermodynamic limit, 
however, the internal energy and entropy described above 
can be calculated equally well from Eqs. ITU)) . We in- 
deed have used these relations and performed the deriva- 
tives with respect to temperature based on finite differ- 
ence approximations using the non-perturbative free en- 
ergies at neighboring temperature as input (see Tab. 1 
of Ref. 0). With this method no perturbative uncer- 
tainties get introduced in our calculations and a test of 
confinement and deconfinement could be given. In fact, 
as we use in our calculations the renormalized free en- 
ergy as input, both quantities, the internal energies and 
entropies, survive the continuum limit. We stress here, 
however, that we also provide non-perturbative renor- 
malization prescriptions for the quark anti-quark inter- 
nal energies and entropies which are independent from a 
renormalization of the quark anti-quark free energies. 

B. Theoretic expectations and renormalization 

Following P| we consider mainly the quark anti-quark 
internal energies and entropies in the color singlet chan- 
nel. At large distances and high temperatures, i.e. 
rT > 1, T sufficiently above T c , as well as at small 
distances and zero temperature, i.e. ^Aqcd 1, the 
singlet free energies, Fi (r, T) , are indeed dominated by 
one gluon exchange [l| ■ Using the perturbative small dis- 
tance relation for the singlet free energy 0, |3(], |H, |3i| , 

Fi(r,T) ~ A^n. , r A QCD «l, (13) 
3 r 

standard thermodynamic relations (Eqs. (JSJ and i|10|)t 

lead to 

C/ 1 (r«l/A QCD ,T) ~ , (14) 

3 r 



1 We assume here the existence of an energy, Vo(r), at T = 
which corresponds to the expectation value of limy_o(C>}. )• 



and 

Si(r« 1/Aqcd,T) ~ 0. (15) 

We thus expect that at small distances the singlet free 
and internal energies are controlled to a large extent by 
energy, i.e. in the limit of small distances both will 
smoothly approach the zero temperature heavy quark 
potential, V(r). At larger distances, however, the quark 
anti-quark free energies are strongly temperature depen- 
dent pi |22i and thus non- vanishing entropy contri- 
butions arise. In this case differences between free and 
internal energies are expected to become important and 
the quark anti-quark free energy will be to a large extent 
controlled by TSi{r, T). 

To clarify the role of the entropy we compare 2 in 
Fig. nja) the short and large distance parts of the sin- 
glet free and internal energies at T ~ 1.3T C and show 
in Fig. [TJb) the corresponding entropy contribution, 
TSi(r,T ~ 1.3T C ). We also indicate in both figures 
the small distance behavior expected from Eqs. 114(1 and 
(|15fl as solid lines, i.e. the line in Fig. [TJa) indicates 
the heavy quark potential from Refs. [ij, l^jj, and in 
Fig. Gib) it indicates the zero level. As both, the in- 
ternal energy and entropy, have been calculated using 
renormalized free energies proper renormalization of both 
observables is already incorporated by construction. It 
can clearly be deduced from Fig. QJa) that the singlet 
free and internal energies smoothly approach V(r) at 
small distances and the free energy thus indeed is dom- 
inated by the energy contribution. In fact, the entropy 
contribution shown in Fig. ^b) is quite small at small 
distances and indicates a vanishing entropy contribution 
in the limit r — > 0. Unfortunately we could not go to 
smaller distances to clearly demonstrate this behavior. 
Moreover, at small distances U±(r,T) and TS\(r,T) suf- 
fer from lattice artifacts which result from small distance 
lattice artifacts in the free energies. At intermediate 
and large distances, however, the free and internal en- 
ergies shown in Fig. UJ[a) deviate from each other and 
TS\ (r, T) indeed plays an important role for the behav- 
ior of F\ (r, T) . At asymptotic large distances the inter- 
nal energies and entropies approach temperature depen- 
dent constant values, i.e. Uoo{T) = lim r _ >00 U\(r,T) and 
Soo (T) = linir^oo Si (r, T) are finite for finite tempera- 
tures. These values are indicated by the arrows in Fig.Q] 
and, in particular, at finite temperature 

E/oo(T) > Foo(T) (16) 

is evident. 

A similar behavior of Uoo(T) and Soo(T) can be de- 
duced from high temperature perturbation theory. To 



2 Details on our lattice simulations, in particular, on the calcula- 
tion of Fi(r, T), are given in Ref. [y . Details on the computation 
of Ui(r,T) and Si(r,T) will be given in Sec. ITTTlandlTVl 
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FIG. 1: (a) The singlet internal energy, Ui(r,T) (filled 
circles), calculated from renormalized singlet free energy, 
Fi(r,T) (open squares), at fixed T ~ 1.3T C calculated in 2- 
flavor lattice QCD and compared them to V(r) (line) [Hl27|. 
(b) The corresponding color singlet quark anti-quark entropy, 
TSi(r,T ~ 1.3T C ), as function of distance calculated from 
renormalized free energies. The arrows point in both figures 
at the temperature dependent values of the free and inter- 
nal energy and entropy at asymptotic large distances, i.e. 
Foo(T) = lim^oo-Fi (r,T), U X (T) = lim r _,oo Ui(r,T) and 
TS^T) = Tlimr-oo Si(r,T). 



be more precise, high temperature perturbation theory 
suggests a cubic leading order dependence 3 of the free 
energy on the coupling 40], i.e. 



Foo(T) ~ --m D (T)a(T) ~ -0(g 3 T) . (17) 



Here a(T) = g 2 {T)/4n and m D (T) denotes the Debye 



3 Here and in what follows we already have anticipated the running 
of the coupling with the expected dominant scale T. Of course, 
the running of the coupling appears only beyond leading order. 



mass which in re-summed leading order is given by 



N, 



m D {T) = (1 + -^) g(T)T. 



1/2 



(18) 



We note that this leading order result is gauge invari- 
ant. In the following we use the renormalization group 
(3- function to evaluate the derivatives of the coupling, i.e. 
for an arbitrary function T = J~(g, T) we use 



T 



dT{g,T) 
dT 



T 



where (3(g) — — /3o<? 3 + 0(g 5 ) in perturbation theory. As- 
suming this behavior the internal energy, Uoo(T), and 
entropy, Soo (T) , are expected to behave like 



U<x{T) ~ Am D (T)a(T) 



/3(g) 
g(T) 



-0(Tg 5 ), (20) 



and 



Soo(T) 



4 mj(T) 
3 T 

-o(g 3 ) ■ 



a(T) 



T 



g(T) 



(21) 



The leading contribution to TS^ (T) is similar to the free 
energy in Eq. I|17|l and thus at high temperatures and 
large distances the free energy is indeed expected to be 
to large extent dominated by the entropy contribution, 
i.e. at leading order Soo(T) ~ -Foo(T)/T. Although the 
entropy itself will vanish logarithmically in the high tem- 
perature limit, i.e. Soq(T — > oo) = 0, the contribution 
TSoo(T) will clearly dominate the differences between 
free and internal energy at high temperatures, 



Uo.iT) - F^T) = TS^iT) 



-0( 5 3 T). (22) 



Thus, the difference between free and internal energy is 
expected to increase continuously with increasing tem- 
perature when approaching the perturbative high tem- 
perature regime. Only in the limit of zero temperature, 
T -> 0, the observable TS , 00 (T) will vanish as S^T) 
is a dimension less quantity which due to string break- 
ing stays finite in QCD. Any qualitative change in the 
observable TS^ (T) as function of temperature between 
both limits, i.e. T — > and T — > oo, is not quite obvious 
and, if present, may signal the phase change from the chi- 
ral symmetry broken phase at low to the deconfincmcnt 
phase at high temperatures. 

We may finally note that although we discuss in the 
following the internal energies and entropies calculated 
from renormalized free energies, it is conceptually quite 
satisfying that both observables could equally well be 
renormalized by matching their short distance parts to 
the heavy quark potential l|14f) and zero l|15|). respec- 
tively. This is indeed evident from Figs.^a, b). As no 
additional divergences get introduced at finite tempera- 
ture also their large distance properties are properly fixed 
in the continuum limit, in particular, also the manifestly 
gauge invariant observables Uoo (T) and (T) . 
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T/T c 




Uoo/T c 


0.79 


5.48 (259) 


9.37 (207) 


0.84 


6.49 (119) 


10.19 (101) 


0.88 


7.78 (168) 


11.30 (148) 


0.93 


12.92 (138) 


15.96 (129) 


0.98 


16.38 (113) 


19.27 (112) 


1.01 


14.83 (172) 


17.71 (173) 


1.04 


12.93 (825) 


15.78 (88) 


1.09 


5.49 (42) 


7.84 (46) 


1.13 


3.95 (41) 


6.14 (46) 


1.19 


2.73 (25) 


4.72 (29) 


1.29 


1.63 (16) 


3.37 (21) 


1.43 


1.25 (11) 


2.85 (16) 


1.57 


1.02 (8) 


2.51 (12) 


1.72 


0.91 (10) 


2.32 (17) 


1.89 


0.87 (7) 


2.26 (14) 


2.99 


0.67 (1) 


2.09 (2) 



TABLE I: The change in internal energy and entropy due to 
the presence of a quark anti-quark pair at infinite separation 
in the QCD heat bath versus temperature. 



III. THE CONFINEMENT DECONFINEMENT 
TRANSITION 

We begin our discussion of the finite temperature en- 
ergies and entropies at asymptotic large distances, r — > 

00. Actually, to avoid any fit we again p| approxi- 
mated the value of the free energy at infinite distance, 
F\(r — > oo, T), by the value of the quark anti-quark 
free energy calculated at the largest possible separation 
on the lattice, i.e. F^T) = F q - q (N a /2,T), and calcu- 
lated separately the internal energy, {/^(T), and entropy, 
Soo(T). Both quantities are obtained from derivatives of 
the color averaged free energy, F qq (r, T) , which is a mani- 
festly gauge invariant observable. Our results for Uoo(T) 
and Sqo(T) are summarized in Tab. [I] It is quite sat- 
isfying that the values obtained for Uao(T) and Soo(T) 
reproduce the free energy given in Tab. 2 of Ref. 

1. e. the quantity Uqo — TSqo matches to the value for 
Foo = — 2Tln|(L}| for all temperatures. 

In parts of our analysis we are also interested in the 
flavor and quark mass dependence of the finite temper- 
ature energies and entropies. We thus again [l] com- 
pare our results (Nf = 2) to results from quenched QCD 
(Nf = 0) 0, EUl and also in parts to a recent study of 
3-flavor QCD |26(. To convert the observables to phys- 
ical units we use T c = 270 MeV in quenched, T c = 200 
MeV in 2-flavor {m v /m p ~ 0.7) and T c = 193 MeV in 
3-flavor {m^/nip ~ 0.4) QCD. It should be obvious, how- 
ever, that a comparison of free and internal energies cru- 
cially depends on the relative normalization of the corre- 
sponding zero temperature heavy quark potentials used 
for renormalization and thus a comparison could be af- 
fected by flavor and/or quark mass dependent (over-all) 



constant contributions. Here, and in what follows, the 
relative normalization of the heavy quark potentials in 
quenched and full QCD is such that there is no constant 
contribution in the Cornell Ansatz for V(r) at large dis- 
tances. We also note that any undetermined constant 
contribution to the heavy quark potential at zero tem- 
perature will add a non-perturbative over-all constant to 
the free and internal energies which would also affect the 
comparison of these observables with perturbation theory 
We stress again, however, that the quark anti-quark 
entropy is unaffected by any undetermined finite renor- 
malization of V(r) at zero temperature, i.e. Soo(T) does 
not dependent on any flavor and/or quark mass depen- 
dent normalization terms that could contribute to V(r) 
at T = 0. 



A. The free energy 

Our results for FodT) are summarized in Fig. as 
function of T/T c and compared to i 7 ' 00 (T) obtained in 
quenched and 3-flavor QCD. While F^T) in quenched 
QCD exhibits a singularity at T c due to the first order 
phase transition and is infinite below T c , it is well-defined 
and finite in full QCD at all temperatures due to string 
breaking below and color screening above the transition. 
In this case Foo {T) is steadily decreasing with increasing 
temperatures in the whole temperature range analyzed 
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FIG. 2: The value of the free energies F X (T) at asymp- 
totic large distances as function of temperature in physi- 
cal units. We compare here our results from 2-flavor QCD 
(T c = 200 MeV, m^/m p ~ 0.7) to the results in pure gauge 
theory (T c = 270 MeV) and 3-flavor QCD (T c = 193 
MeV, rn^/nip ~ 0.4) |2ql . The dashed horizontal lines show 
the expected energy V(rbreaking) — 1000 ~ 1200 fm using 
t breaking — 1.2 — 1.4 fm from 2-flavor lattice studies at T = 
and quark mass m^jm p ~ 0.7 43]. In 2-flavor lattice studies 
at T = and lower quark mass, m^-/m p ~ 0.4, the string is 
expected to break at smaller (lower) distances (energies) [4411 . 
The thick line around T c is explained in Sec. 1111 51 
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FIG. 3: The contribution TS'oo(T) appearing in the free en- 
ergy, F X (T) = Uoo (T) — TSao (T) , calculated in 2-flavor QCD 
as function of T/T c . We compare our results from 2-flavor 
QCD to the leading order part of Eq. (or Eq. l(2*2*) 1 us- 
ing the 2-loop formula for the coupling 1231 with Xc/Ajjg = 
0.77(15) [Uli^ and scales n = vr, ...,4tt (dashed lines). We 
also show results from 3-flavor QCD for T > 1.1T C |2ql and 
quenched QCD HEH^. 



by us. A discussion of (T), in particular for T < T c , 
has already been given 0, |23| ■ We may add here that 
the slope of Foo(T) as function of temperature indeed 
turns out to be maximal in the vicinity of the transition. 
Although flavor and/or quark mass dependences of the 
observable F^iT) can clearly be seen when comparing 
the data for 2- and 3-flavor QCD, at temperatures about 
1.2T C < T < 2T C no or only little differences between 
quenched and full QCD could be identified. Flavor de- 
pendences are, however, to be expected when reaching 
temperatures in the perturbative regime. 



B. The entropy 

The entropy contribution, T5 00 (T), obtained in 2- 
flavor QCD is shown in Fig. [3] as function of T/T c and 
is again compared to results from quenched and 3-flavor 
QCD. We indeed find T5 , 00 (T) at all temperatures 
analyzed here. Moreover, the data for TS' 00 (T) at low 
temperatures also suggest a vanishing contribution in the 
zero temperature limit, T — > 0. Unfortunately we could 
not go to smaller temperatures to clearly demonstrate 
this behavior. On the other hand, in the high tempera- 
ture phase, i.e. at T > 2T C , we indeed find a tendency 
for an increase of TSqo (T) with temperature as expected 
from i|22|) . Actually, this small increase is consistent with 
the rise given in (|2ip. To demonstrate this we also com- 
pared TS , 00 (T) in 2-flavor QCD to Eq. using the 



perturbative 2-loop coupling, i.e 



92-ioop(T) = 2/3 ln 



A 



MS 



MS 



(23) 



with 



0o 

ft 



l 



16tt 2 
1 



11 - 



2N, 



(16tt 2 ) 5 



102 



38N, 



assumin g va nis hing quark masses. We used T c /Aj^g = 
0.77(15) TUEIIll and the ambiguity in fixing the scale 
in perturbation theory, (i — n, ..., Ait. This estimate is 
shown within the dashed lines and qualitatively agrees 
with the lattice data for T > 2T C . However, to clearly 
establish the perturbative increase of T5 , 00 (T) with in- 
creasing temperature will require the analysis of signif- 
icantly higher temperatures. We also note that within 
the statistical accuracy of the data we find for T > 2T C 
the tendency, 



S 



Nj=0 



(T) < S^~ 2 (T) < S%- 3 {T) 



(24) 



It appears indeed quite reasonable that introducing ad- 
ditional flavor degrees of freedom may enhance the finite 
temperature quark anti-quark entropy in the deconfmed 
phase. It is, however, quite difficult to separate clearly 
the different effects from flavor and quark mass depen- 
dence in full QCD. In particular, at temperatures below 
T c the tendency given in \2ty may change as can be seen 
from the temperature dependence of FociT) shown in 
Fig.El 

In contrast to the small temperature dependence of 
T5 , 00 (T) at low and high temperatures, T5oo(T) shows 
qualitatively and quantitatively significant differences at 
temperatures in the vicinity of the transition. In fact, 
TSooiT) obtained in 2-flavor QCD exhibits a sharp peak 
at T c . This behavior signals the high temperature phase 
transition/crossover in QCD. As the peak is so sharp we 
may introduce TJ (T u ) defined as the lower (upper) tem- 
perature at which Soa(T) approaches about half of the 
peak value, i.e. Soo (?/,„) = S QO (T c )/2. We find for 2- 
flavor QCD T x about 0.89T C and T u about 1.07T C using 
Soa(T c ) — 16.5. This temperature range is shown at the 
bottom of Fig. [21 as thick line. A similar behavior is 
also apparent in 3-flavor QCD. This indicates that the 
crossover from the low to the high temperature phase in 
QCD takes place in a small temperature range around T c . 
We stress, however, that the temperature dependence of 
Foo(T) also at temperatures above 1.07T C is still to a 
large extent dominated by non-perturbative effects. 



C. The internal energy 

The internal energy, [/^(T), in 2-flavor QCD is shown 
in Fig. 0| as function of temperature and is compared 
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FIG. 4: The internal energy U X (T) versus T/T c calculated 
in 2-flavor QCD. The corresponding free energy, Fao(T), cal- 
culated in 2-flavor QCD is also shown as solid line. We again 
indicate in this figure the energy at which string breaking is 
expected to take place at T = 0, V^rbreaking) — 1000 — 1200 
MeV (dashed lines), using rbrcaking = 1-2 — 1.4 fm [JJ. 



to the corresponding free energy, F^T) (solid line), al- 
ready shown in Fig. □ We indeed find U^T) > F^T) 
at all temperatures analyzed here. It can clearly be seen 
that the temperature dependence of (T) and Uoc (T) is 
qualitatively and quantitatively different. While the free 
energy steadily decreases with increasing temperatures 
the internal energy exhibits a pronounced peak. Again 
this peak is sharply localized at the (pseudo-) critical 
temperature. At the temperatures analyzed by us the 
internal energy below T c is rapidly increasing with in- 
creasing temperatures. Again we indicate by dotted lines 
the plateau value of the heavy quark potential at zero 
temperature, F(rbreaking) — 1000 — 1200 MeV, using 
^breaking — 1-2 — 1.4 fm ^ . A comparison of Uoc(T) 
with this value shows again that most of the temperature 
dependence of Uoo(T) is sharply localized at tempera- 
tures in the vicinity of the transition. A qualitatively sim- 
ilar behavior is also apparent in 3-flavor QCD 26]. Com- 
paring the available data some flavor or quark mass de- 
pendence, uH i=2 (T c ) ~ 4000 MeV > u" f=3 (T c ) ~ 3000 
MeV, can be observed and also the value J7 00 (T C + ) in 
quenched QCD ^l| is of similar magnitude. As noted 
in Sec. Ill Bl however, the values for Uoo(T) may depend 
on the relative normalization of V(r) at T = used for 
renormalization. 

At temperatures above T c , Uoo(T) rapidly drops while 
at higher temperatures, i.e. T > 1.3T C , the tempera- 
ture dependence of Uoo(T) turns out to be much weaker 
than in the vicinity of phase transition. However, con- 
tact with the perturbative relation, Eq. is not ex- 
pected at those temperatures as Uoo(T) is still positive. 
In fact, when approaching the perturbative high temper- 
ature regime also Uoq (T) is expected to exhibit a change 
in sign and will slowly diverge with respect to l|20|) . 



IV. FINITE TEMPERATURE ENERGIES AND 
POTENTIAL MODELS 

The analysis of bound state problems has bee n q uite 
successful in terms of potential theory at T = |48t I4H 
50] . For the discussion of quarkonium suppression pat- 
terns at finite temperature one also often resorts to po- 
tential models [fll fllL HH 111 Il7j . Of course, the strong 
interaction remains unaffected by temperature and the 
modeling of thermal modifications of heavy quark bound 
states requires the definition of an effective potential, 
V e g(r, T) 0, which can be given only phenomenologi- 
cally, for instance, by using the modifications of the free 
and internal energies. It is thus important to understand 
the binding properties of the different finite temperature 
energies. For this purpose we also calculated the quark 
anti-quark internal energies for the temperatures given 
in Tab. U at several finite distances. Parts of our results 
for Ui(r,T) are summarized in Fig. at temperatures 
below (a) and above (b) the transition. Similar results 
have been obtained for internal energies in the averaged 
and octet channels. 

To gain some insight into the consequences these en- 
ergies have for quarkonium dissociation we compare the 
asymptotic (r — » 00) energies with the potential energies 
at a distance corresponding to the size of some quarko- 
nium states, 

AEi(T = 0) = V(oo) - V(n) , (25) 

where the radii r, (A — J/ip, Xc---) are listed in the 
first row of Tab. [H] At zero temperature V^(oo) is 
taken to be twice the energy needed to create the low- 
est heavy-light meson. For our purpose we consider 
V(oo) = ^(rbrcaking) ^ 1000 MeV where rebreaking is 
the distance at which the string is expected to break 
at zero temperature, rb rea king 1.2 fm 43]. This en- 
ergy is shown in Fig. as horizontal line. The resulting 
energy for J/ip is also shown. The energies for some 
charmonium and bottomonium states are summarized in 
Tab. In] and compared to the mass difference obtained 
from 2Md : b ~ mi where Md,b denotes the D- and B- 
meson masses and rrii the masses of the different quarko- 
nium states jU- Of course, the wave functions for the 
different quarkonium states will also reach out to larger 
distances [5(| and thus our estimate for the different en- 
ergy levels Ei (T = 0) can only be taken as indicative for 
the relevant energies. Potential model analysis, using for 
instance the Schrodinger equation, will do better in this 
respect. 

Similarly we can estimate the temperature dependence 
of the energy levels for the different quarkonium states 
from Ei(T) = V c s(ri,T). Again these energy levels will 
only characterize the relevant energies and the sizes of 
these states may also become temperature dependent. 
At finite temperature, however, the values for these levels 
are expected to depend crucially also on the specific mod- 
eling of the effective potential, V e g(r,T). This is obvious 
from the different energy levels for the J/ij) shown in 
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state 


J/i> 


Xc 


ip> 


T 


Xb 


T/ 


n [fm] 


0.45 


0.70 


0.87 


0.23 


0.41 


0.51 


2M d ,b - rm [MeV] 


631 


313 


42 


1098 


698 


535 


AEi(T = 0) [MeV] 


830 


560 


400 


1140 


870 


750 


AEi(T c ) [MeV] [using Fi] 


321 


124 


62 


665 


334 


265 


AEi(T c ) [MeV] [using Ui] 


3036 


1618 


1010 


3721 


3085 


2615 


AEi(T c ) [MeV] [using F qq ] 


97 


27 


12 


607 


87 


55 


AEi(T c ) [MeV] [using Uq q ] 


1355 


747 


504 


3732 


1392 


1136 


Id™ /r c [potential model using Fx] 


1.1 


0.74 


0.1 - 0.2 


2.31 


1.13 


1.1 


T™/T c [potential model using Ui] 


~ 2 


~ 1.1 


~ 1.1 


~4.5 


~ 2 


~ 2 


from lattice spectral functions 


1.5 - 3 


< 1.1 







TABLE II: A summary of the different estimates for the relevant temperatures, Tdis, at which quarkonium dissociation is 
expected to become important. In the first row we give the mean squared charge radii, ri, and the energies, AEi(T = 0) = 
V(oo) — V(ri) at zero temperature with V(oo) ~ 1000 MeV. We also list the mass gap 2Mb,b — rm using D- and B-meson 
masses M D = 1864 MeV and Mb = 5279 MeV, respectively (EH- The second row summarizes results for the break up energies 
at T c estimated from AEi(T c ) = V c s(ri,T c ) — V c g(oo,T c ) using different definitions for V e «(r,T). The errors on these values are 
typically about 20%. Note here that the thermal energy E th {T) = 3T is about 600 MeV at T c . The 3 rd and 4 th row contain 
different values for the relevant dissociation temperatures, i.e. estimates for the onset of temperature effects in quarkonium 
states, T nset 0-0 the predicted dissociation temperatures from potential model calculations (TJis^) fil H3. Il4|. and lattice 
studies of charmonium spectral functions, Tj^ |52ll53||. 



Fig. |5]which we obtained by using as V c s(r, T) the singlet 
free energy (lower dashed line) and singlet internal energy 
(upper dashed line). Due to the steeper rise of the inter- 
nal energy compared to free energy Ej/^(T) is enhanced 
compared to the energy level obtained from the inter- 
nal energy. It is interesting to note here that Ej/^(T — 
1.3T C ) deduced from the internal energy is even larger 
than at zero temperature while in terms of the free energy 
it is smaller than Ej/^(T = 0). For the characterization 
of the relevant energies needed to dissociate the bound 
state we again consider A£!j(T) = V^g(oo, T) — V c s(ri,T), 
which depends on the definition of V c g(r, T) . This is evi- 
dent from Fig. [7| where we show the temperature depen- 
dence of AEj/^T) estimated from V cS (r,T) = U x {r,T) 
(filled circles) and V e g(r, T) = F\{r,T) (open circles). 
While AEj/^(T) is continuously decreasing with increas- 
ing temperatures when using F\ (r, T) as effective poten- 
tial, its binding pattern appears quite different from what 
one obtains by using V e s(r, T) = U±(r,T). Actually, in 
the latter case the binding pattern exhibits a maximum in 
the vicinity of the transition while it rapidly drops above 
T c . Similar results have been obtained also for other char- 
monium and bottomonium states and are summarized in 
Tab.HU 

Of course, the model dependences of the dissociation 
energies at finite temperature also affect the analysis of 
suppression patt erns and corresponding dissociation tem- 
peratures H1H0. Actuall y> usin S Ves(r,T) = Ux{r : T) 
suggests that suppression of J/tp may occur only at 
temperatures close but above the transition while from 
V e g(r, T) = F\ (r, T) one finds that J/ tfj dissolves already 
at temperatures below the crossover. Similar model de- 
pendences enter also the analysis of excited quarkonium 
states and the corresponding estimates for the dissocia- 



tion temperatures are summarized in Tab. ITT1 using four 
different definitions for V c g(r,T), i.e. we used the sin- 
glet free and internal energies (F\(r, T), U\(r, T)) as well 
as the finite temperature energies in the color averaged 
channel (F qq (r,T), U qq (r,T)). 

V. SUMMARY AND CONCLUSIONS 

Following we introduced and analyzed the change 
in internal energy and entropy due to the presence of 
a static quark anti-quark pair in a QCD heat bath. 
Both observables are introduced as intensive observables 
as appropriate derivatives of the renormalized free en- 
ergy. Similar to the singlet quark anti-quark free ener- 
gies pi 0, ^2 El also the singlet internal energies become 
temperature independent in the limit of small distances 
and are controlled by the zero temperature running cou- 
pling. 

We analyzed qualitative and quantitative differences 
that appear when changing from free energies to inter- 
nal energies as observable that defines an effective po- 
tential that can be used in model calculations. We dis- 
cussed the important role of the entropy contribution 
at finite temperature. At short distances, at interme- 
diate and at large distances the entropy contribution 
is non-zero and shows non-trivial r-dependences. We 
find positive entropy contributions, S\(r,T) > 0, and 
thus Ui(r,T) > Fi(r,T). Similar to the free energies, 
also the large distance properties of the internal ener- 
gies and entropies are controlled by string breaking be- 
low and color screening above deconfinement and both 
approach temperature dependent constant values, which 
define Uqo (T) and Soo (T) , at asymptotic large distances. 
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observables which can indeed be seen after renormaliza- 
tion in the temperature dependence of the finite temper- 
ature energies, entropies and the renormalized Polyakov 
loop (see also PHEE3)- In contrast to quenched QCD, 
in full QCD the phase change is a crossover and we con- 
sequently do not see any singularities in the finite tem- 
perature energies and entropies nor in the temperature 
dependence of the renormalized Polyakov loop [lj . 

We also investigated the temperature dependence of 
quarkonium binding in the vicinity of the transition. For 
this purpose we used the finite temperature energies to 
define appropriate effective potentials, T4ff(r, T), as one 
would do in potential models. As effective potentials 
we used the singlet and averaged free and internal ener- 
gies, i.e. we defined V e s(r, T) through Fi(r, T), Fq q (r, T), 
U\{r,T) and Uq q (r,T), and estimated the binding ener- 
gies at temperatures below and above the transition. In 
all cases the binding energies of the quarkonium states in- 
deed become weaker with increasing temperatures above 
the transition and this may lead to dissociation of parts 
of these states at temperatures close but above T c . Our 
analysis, however, shows strong dependencies of the bind- 
ing energies and dissociation temperatures on the spe- 
cific modeling of V e s(r,T). To some extent these model 
dependencies enter from quite general grounds when ex- 
changing the definition of V c g (r, T) from free energies into 
internal energies. When using free energies the binding 
energies continuously and rapidly decrease when crossing 
the transition and most of the quarkonium bound states 



FIG. 5: The color singlet quark anti-quark internal energies, 
Ui(r,T), at several temperatures below (a) and above (b) the 
phase transition obtained in 2-flavor lattice QCD. In (a) we 
also show as horizontal lines the asymptotic values given in 
Tab. E which are approached at large distances and indicate 
the flattening of Ui(r,T). The solid lines represent in each 
figure the T — heavy quark potential, V(r) [ILl27||. 



Actually, the difference between free and internal ener- 
gies at high temperatures, TS , 00 (T) ~ 4m D (T)a(T)/3, 
is supposed to increase with increasing temperatures. In 
particular, Uoo(T) and Soq(T), are, similar to Foo(T) ||, 
again introduced as manifest gauge invariant observables 
and clearly signal the QCD plasma transition. In fact, 
while the plateau values which are approached by the free 
energies, FooCT), are rapidly decreasing in the vicinity of 
the transition [2^ , the values approached by the inter- 
nal energies, Uoo(T), and entropies, Soo(T), show both a 
sharp peak at the (pseudo-) critical temperature. Similar 
results are also obtained in quenched and 3-flavor QCD 
[E IE HE EH- However, qualitative differences become 
quite transparent in the vicinity and below the transition 
when comparing these observables obtained in quenched 
and full QCD. In quenched QCD the first order phase 
transition is related to singularities in thermodynamic 



V eff (r,T) [MeV] 




0.5 1 1.5 



FIG. 6: Different effective potentials in the color singlet chan- 
nel, Vrff(r,T) = Fi(r,T) (open symbols) and V e g(r,T) = 
Ui(r,T) (filled symbols), at fixed T ~ 1.3T C as function of 
distance and the heavy quark potential, V[r). We also com- 
pare V e g(r, T) to the J/ip energy level at T = 0, Ej/^(T = 
0) = V(r = rj/ij,) ~ 270 MeV (horizontal solid line). The 
horizontal dashed lines correspond to the J/ip energy levels 
defined on Ej/^(Fi) = Fi(rj/i/j,T) (lower dashed line) and 
Ej/^(Ui) = Ui(rj/^,T) (upper dashed line), rj/^ ~ 0.45 fm 
is fixed through the mean squared charge radius expected at 
T = given in Tab.|TT| 
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increasing binding energies below the transition which 
exhibit a peak at T c . This may imply that all quarko- 
nium states analyzed here are still bound at the transi- 
tion. The binding energies of the different states rapidly 
decrease above T c leading again to quarkonium dissoci- 
ation, however, the temperatures which are relevant for 
dissociation are shifted to larger temperatures than those 
deduced from free energies. The recent potential model 
calculations 0, using the properties of U\ (r, T) at 
temperatures above the transition as well as the lattice 
analysis of quarkonium bound states in quenched QCD 
|52ll53j| support our findings. 
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FIG. 7: The energy difference V cS (oo,T) - V cB (rj/^,T) 
as function of temperature using the differ- 
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